THz Generation and Detection on Dirac Fermions in Topological Insulators 
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This study shows that a terahertz (THz) wave can be generated from the (001) surface of cleaved 
Bi2Se3 and Cu-doped Bi2Se3 single crystals using 800 nm femtosecond pulses. The generated THz 
power is strongly dependent on the carrier concentration of the crystals. An examination of the 
dependence reveals the two-channel free carrier absorption to which Dirac fermions are indispensable. 
Dirac fermions in Bi2Se3 are significantly better absorbers of THz radiation than bulk carriers at 
room temperature. Moreover, the characteristics of THz emission confirm the existence of a recently 
proposed surface phonon branch that is normalized by Dirac fermions. 

PACS numbers: 71.27.+a, 78.47.D-, 78.68.+m 
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Three-dimensional topological insulators (TIs) are 
characterized by a narrow band gap in the bulk and a 
Dirac cone-like conducting surface state [H-Q. The sur- 
face state is a new state of quantum matter caused by 
the strong spin-orbit interaction and protected by time- 
reversal symmetry. The special properties of TIs have 
applications in spintronics and quantum computations. 
Certain TIs with a small band gap are especially use- 
ful for terahertz (THz) optoelectronics. One of the key 
issues about TIs has been the identification of the gap- 
less surface electronic states (Dirac fermions) and the 
characterizations of their fundamental properties. Angle- 
resolved photoemission spectroscopy (ARPES) [l|, |j-|6[ 
and scanning tunneling microscopy (STM) jTHlOf have 
successfully confirmed the existence of Dirac fermions 
in Bii.xScx, Bi2Se3, and Bi2Te3. Regarding the trans- 
port measurements, a metallic channel associated with 
the protected surface state has been detected by cither 
controlling the gate voltage in TI devices with a suffi- 
ciently low bulk carrier density or by using very thin TI 
films [ll-14|. However, these experiments used partic- 
ularly specialized instruments, and their procedures are 
excessively complex for quick and routine characteriza- 
tions of Dirac fermions in TIs. Hsieh et al. [l^ showed an 
alternative approach using second harmonic generation 
(SHG) in arsenic- doped Bi2Se3 single crystals associated 
with Dirac fermions, which showed a new venue for exam- 
ining Dirac fermions by contact-free optical techniques. 
However, SHG is highly sensitive to the surface quality 
of samples and doping. THz waves may in principle be 
an ideal tool for distinguishing Dirac fermions from bulk 
carriers because they are not sensitive to the surface qual- 



ity of samples with long wavelengths. Furthermore, THz 
wave have a photon energy (approximately 4 meV) that is 
significantly lower than the bulk gap (approximately 300 
meV) of TIs; thus, THz radiation would allow specific 
characterizations within the Dirac cone. Aguilar et al. 
[lij recently showed the THz responses of Dirac fermions 
in Bi2Se3 thin films. However, this type of THz exper- 
iments can only be applied to thin TIs of several tens 
of quintuple layers. This study shows a THz generation 
from pure Bi2Se3 and Cu-doped Bi2Se3 single crystals by 
pumping with femtosecond laser pulses. Dirac fermions 
were identified to have an indispensable role on the inten- 
sity of THz emission. Free carrier absorption is a crucial 
mechanism to the optoelectronic devices of TIs and was 
revealed from the dependence of generated THz power 
on the carrier concentration 
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Moreover, the detailed 
characteristics of THz generation verified a newly pro- 
posed phonon branch normalized by Dirac fermions jl8| . 

Single crystals of pure Bi2Se3 and Cu-doped Bi2Se3 
were grown using either the Bridgeman, Melt growth, or 
CVT methods [19[. Single crystals of CuxBi2Se3 were ob- 
tained using a slow-cooling method from 850 to 650 °C 
at a rate of 2 °C/h and quenching in cold water. Scotch 
tape was used to cleave the (001) surface of the Bi2Se3 
crystals to ensure a flat and bright surface for optical 
measurements. The carrier concentrations of the samples 
listed on Table I were obtained using the Hall measure- 
ments. The mobility was measured using the four-probe 
method. A reflection-type THz generation scheme was 
used to generate a THz wave on TIs, as shown in Fig. 
[IJa) and the inset of Fig. [IJb). An 800 nm Ti:sapphire 
laser (FemtoLasers, Inc.) beam with a repetition rate of 




TABLE I. Carrier concentration and the THz peak amplitude 
for samples grown by diflFerent methods Bridgman, CVT, 
Melt growth). All samples are n-type. 



Carrier 

Code Compounds concentration 



(-10^ 



THz peak amplitude 
(arb. units.) 



"#1 
^#8 



BiaSeg 

BiaSeg 

BiaSea 

BiaSea 
Cuo.oaBiaSes 
Cuo.osBiaSes 
Cuo.iBiaSes 
Cuo.i25Bi2Se3 



75.5±13.6 
34.6±4.37 
31. Oil. 61 
15.6±10.3 
3.66±0.16 
4.23±1.17 
1.96±0.86 
1.17±0.56 



0.72±0.62 
4.49±0.67 
3.05±0.51 
6.43±0.43 
31.5±0.44 
32.3±0.51 
22.8±0.17 
30.3±0.22 



5.2 MHz and a pulse duration of 50 fs was incident at 9 
= 45° (to the surface normal) and focused on the surface 
of the samples with a diameter of 43 /zm. The pumping 
fluence was tuned by varying the laser output power (the 
typical value for this study was 0.37 mJ/cm^). Follow- 
ing femtosecond pulse pumping, the generated THz wave 
was collected using a pair of off-axis parabolic mirrors 
and focused on a 1-mm-thick ZnTe crystal to allow its 
detection with electro-optical (EO) sampling [2^. The 
entire generation and detection systems were sealed in a 
nitrogen-filled plastic box to reduce the humidity to < 
6.0%. All optical measurements were performed at room 
temperature. 

Fig. lUb) shows the typical THz waveform generated 
from pure Bi2Se3 and Cu-doped Bi2Se3 single crystals 
with a reflection-type setup (THz radiation cannot be de- 
tected after TIs). The amplitude of the THz wave gen- 
erated from pure Bi2Se3 single crystals is significantly 
smaller than that from a Cuo.o2Bi2Se3 single crystal. 
Certain Bi2Sc3 crystals such as sample #1 produce nearly 
zero amplitude (below the S/N ratio in the detection sys- 
tem). The THz generation intensity is strongly depen- 
dent on carrier concentration and doping. In general, the 
THz waveform is composed of a large single pulse and a 
damped oscillation, which is due to the interference be- 
tween THz electric fields from different positions of the 
crystal. These electric fields are caused by the mismatch 
between the THz phase velocity and the group velocity of 
the optical pumping pulse [2lj . Time-domain THz wave- 
forms (Fig. [Ijb)) can be converted to frequency domain 
spectra (Fig. [2|) by using fast Fourier transform (FFT). 
The central frequency for Cuq. 02612863 and bandwidth is 
approximately 1.2 THz and 1.6 THz, respectively. The 
intensity of pure Bi2Se3 THz spectra is relatively small, 
corresponding to small THz signals in the time domain. 

THz signals were measured at various azimuth angles tp 
along the surface normal to understand the THz genera- 
tion mechanism in TIs (inset of Fig. [TJb)). The data from 
inset of Fig. [3]clearly show that the THz peak amplitude 
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FIG. 1. (color online) (a) Schematic illustration of THz gen- 
eration and FCA inside topological insulators. The dots in- 
dicate the free carriers, (b) THz waveform generated from 
various BiaSea and Cuo.oaBiaSea single crystals, and an InAs 
wafer. Inset: schematic illustration of THz generation and 
the detection scheme. 



is virtually independent of (p and at odds with the opti- 
cal rectification simulation curve with six-fold symmetry 
[2^ . Therefore, optical rectification is not the dominant 
mechanism for THz generation in TIs, and the nonlinear 
effect does not mainly contribute to THz generation in 
TIs. 

The band gap of 0.3 eV in Bi2Se3 is significantly 
smaller than the pumping photon energy of 1.55 eV. Free 
carriers are generated when the femtosecond laser illu- 
minates Bi2Se3 or Cu-dopcd Bi2Se3 crystals. These ex- 
cited carriers are located inside the bulk within 100 nm 
[isl ]. When the electric field is built inside the crystals, 
the excited carriers in the bulk are driven and form the 
currents. Two types of built-in electric fields are nor- 
mally present in semiconductors. The surface depletion 
field results from the bending of the conduction band on 
the semiconductor surface as in GaAs [l^]. The photo- 
Dember field is caused by the inhomogeneous distribution 
of holes and electrons [25|, which is usually present in 
narrow bandgap semiconductors, such as InAs and InSb 
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FIG. 2. (color online) Fourier power spectra in the frequency 
domain converted from the time-domain THz signals in Fig. 
\V[h) using fast Fourier transform for Bi2Se3 and Cuo.o2Bi2Se3 
single crystals, and an InAs wafer. Inset: the THz peak am- 
plitude of a Cuo.o2Bi2Se3 single crystal as a function of the 
pumping fluence. 




Frequency (THz) 

FIG. 3. (color online) Polarization-dependent Fourier power 
spectra of a Cuo.o8Bi2Se3 single crystal. Pthz (Sthz): The 
electric field of the p-(s-)polarized THz radiation is parallel 
(perpendicular) to the plane of incidence (see the inset of Fig. 
[TJb)). The solid squares are the surface topological phonon 
dispersion data taken from Ref. [31 • Inset: THz peak am- 
plitude for a Cuo.o2Bi2Se3 single crystal as a function of the 
azimuth angles (fi along the [001] direction, as shown in the 
inset of Fig. [Ub). 



with a bandgaps of 0.35 cV and 0.17 cV, respectively. 
The bandgap is approximately 0.3 eV for Bi2Se3 single 
crystals, which is close to that of InAs (reference sample 
in Fig. [T]). Additionally, the intrinsic charge inhomogene- 
ity in the vicinity of the surface and, result, band- 
bending effects were reported [26]. Consequently, both 
photo-Dember and surface depletion effects are possible 
mechanism for THz generation in Bi2Se3 and Cu-doped 
Bi2Se3 single crystals. The currents formed by the ex- 
cited carriers are suppressed within several picoseconds, 
due to carrier scattering with impurities (e.g., Se vacan- 
cies) or with the layer boundary. The transient current 
further generates THz radiation by ExHzit) oc dJ{t)/dt. 
Higher pumping fluences generate more free carriers, 
leading to larger changes of the transient current; thus 
, a stronger THz radiation should be generated. The 
THz peak amplitude increases linearly with the pumping 
fluences (inset of Fig. ^ 

The excited carriers in the bulk can diffuse either along 
the [001] direction or on the (001) plane to form two 
types of currents. In principle, the diffusion along the 
[001] direction is suppressed by the layer boundary or im- 
purity scattering to generate p-polarizcd THz radiation, 
and the diffusion on the (001) plane is suppressed mainly 
by impurity scattering to generate s-polarized THz ra- 
diation. The addition of a single wire-grid polarizer be- 
tween the sample and an EO detection system allows the 
p- and s-polarized THz radiation from TIs (e.g., from 
the Cuo.o8Bi2Se3 single crystal) to be distinguished (Fig. 
|3]). Intriguingly, the difference between p-polarized and 
s-polarized THz radiation is not only in the central fre- 
quency, but also in the shape of the spectra. Subtraction 
of the s-polarizcd THz spectrum from the p-polarizcd 



THz spectrum results in the gray area between 0.53 THz 
and 1.72 THz (Fig. [Sj, indicating additional absorp- 
tion for the s-polarized THz wave. Zhu ct al. [l^ re- 
cently measured the surface phonon dispersion in Bi2Se3 
crystals using the coherent helium beam surface scatter- 
ing technique. They discovered a low-energy isotropic 
convex-dispersive surface phonon branch normalized by 
Dirac fermions and with a frequency range from 1.8 to 
0.74 THz. This frequency range is consistent with that of 
the gray area in Fig. [3l Therefore, the missing power of 
an electromagnetic wave with an electrical field parallel 
to the surface manifests the additional effect due to the 
Dirac-fermion-normalized surface phonons. 

Table I shows that the carrier concentration decreases 
by more than one order of magnitude when Cu is 
doped into Bi2Se3 crystals (n-type, caused by Se va- 
cancies) because Cu replaces Bi in the lattice 28-30|. 
The carrier concentrations from 15.6±10.3xl0^^ cm"^ to 
75.5±13.6x 10^^ cm"^ were observed for the pure Bi2Se3 
crystals, potentially because of the different growing con- 
ditions and methods. The THz signals from pure Bi2Se3 
crystals are generally relatively small. Conversely, the 
Cuo.o2Bi2Se3 crystal with a lower carrier concentration 
produces stronger THz emission than pure Bi2Se3 crys- 
tals. The spectral weight of the FFT spectra in Fig. 
[2] was plotted as a function of the carrier concentration 
to further quantify these results. Fig. U clearly shows 
that the THz output power (i.e., the spectral weight of 
the FFT spectrum) increases with the decreasing carrier 
concentration. 

The THz wave generated inside the bulk likely suffers 
free carrier absorption (FCA) during propagation. This 
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effect can generally be reduced by suppressing the carrier 
concentration to increase the output intensity of the THz, 
as demonstrated in III-V semiconductors such as InAs 



31| . The dependence of the THz intensity on the carrier 
concentration can be phenomenologically described using 
the Beer-Lambert Law: 



Ithz — -?0,THze 



— crln 



(1) 



where /thz and /q.thz are the transmitted THz intensi- 
ties outside and inside the samples, respectively, as shown 
in Fig. [IJa), I is the path length, n is the carrier concen- 
tration, and a is the absorption cross-section. The exper- 
imental data in Fig. |4]can be fitted as the black dashed 
line by Eg. ([T]), and is qualitatively similar to the case 
of InAs [31[. However, a closer inspection reveals that 
the fit by Eq. ([1} leads to a larger deviation from the 
data in the low concentration regime. Consequently, the 
experimental data in Fig. |4] cannot be explained solely 
by FCA due to the bulk carriers. 

It is known that Dirac fermions exist in Bi2Se3 crys- 
tals P . The Dirac band contribution to FCA on the THz 
emission must be considered. Fig. [TJa) shows that the 
THz wave can be absorbed by the electrons on the Dirac 
cone. The effective two-channel FCA including the con- 
tribution from both the bulk carriers and Dirac fermions 
is written as a modified Beer-Lambert equation: 



(2) 



where erf, is the absorption cross-section of bulk, lb is the 
path length of the bulk, n;, is the carrier concentration 
in bulk, (7s is the absorption cross-section of the surface 
state. Is is the path length of the surface state, and Us 
is the carrier concentration in the surface state. The 
empirical relation between rib and can be described 
by nsXlO-^3=o.51+2.10[l-exp^ri6x/bX 10-13/20.48)] as 
shown in the inset of Fig. 2] [32[. The observed spec- 
tral weight of THz in Fig. 2] is significantly better fit by 
Eq. ([5]) (the red solid line) with Ig = 2 nm (thickness 
of the surface state) [s^ and lb = 23.5 nm (THz emit- 
ted from the bulk within 23.5 nm; i.e., the penetration 
depth of 800 nm pumping light) [34| . The acquired fitting 
parameters are <7b =2.59x10"^^ cm^ and CTs =1.46x10"^'^ 
cm^, and the ratio of Osjob = 5.64. Therefore, this study 
concludes that FCA caused by Dirac fermions is more ef- 
ficient than that caused by bulk carriers (at least at room 
temperature). 

FCA can be explained in the context of the Drude 
model. The FCA cross-section 
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(3) 



where e is the electron charge, [l is the mobility, w is 
the angular frequency of the THz wave, t is the scatter- 
ing time of the carriers, n is the refractive index of the 
material, c is the seed of light, and eg is the permittiv- 
ity constant. Equation ^ indicates that the observed 
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FIG. 4. (color online) Spectral weight of the fast Fourier 
transform (FFT) spectra vs. the carrier concentration. The 
black dashed line represents the fit by Eq. U]). The red 
solid line represents the fit by Eq. ((2]). Inset: the car- 
rier concentration ris of the Dirac fermions vs. nj, of the 
bulk carriers. The dotted line represents the empirical fit 
of n,xl0"^=*=0.51+2.10[l-exp(-n6x/6Xl0"^V20.48)]. 



cross-sectional ratio is as follows: 



Ob 



Ms (1 



M6 (1 



(4) 



where s denotes the physical properties of the Dirac 
fermions of the surface state and h denotes the physi- 
cal properties of the bulk carriers. In this study, /if, « 
1500 cm^/V s at room temperature for the typical Bi2Se3 
single crystals. The mobility can be expressed as \ib = 
eT/m, where m is the effective mass of the carriers. With 
the effective mass of the bulk carrier mb=0-l'imo (35l.[36t. 
where mo is the electron bare mass, Tb ~0.12 ps at room 
temperature. According to Ref 36 and the references 
therein, ^ Tf, at low temperatures (low T); so is it as- 
sumed at room temperature. The above analysis results 
in ujTb ~ 1.2 and loTs 1. Usually, is smaller than fib 
at low T in Bi2Se3 single crystals of approximately lOO 
fim in thickness [ill \3aL [36j . However, the experimen- 
tal data for fj,s at room temperature have been relatively 
rare. Equation Q indicates that the ratio (Js/crb ~ 5.64 
may lead to comparable values for fig and /Xb at room 
temperature. Generally, the effective mass varies weakly 
with T. Conversely, Tb decreases rapidly with increasing 
T because of electron-phonon scattering. The scattering 
of Dirac fermions may be more dominated by the im- 
purity scattering than that of the bulk carriers; thus, 
is less temperature dependent and does not decrease as 
rapidly with increasing T compared to Tb- Therefore, the 
values of fis and fib may be comparable at room temper- 
ature. These results strongly suggest that a THz wave 
generated inside the Bi2Se3 crystals can be used not only 
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to identify but also to effectively examine fundamental 
properties of Dirac fermions. 

In summary, THz radiation can be generated from 
Bi2Se3 and Cu-doped Bi2Se3 single crystals. Dirac 
fermions of the surface state are indispensable to explain- 
ing the strong dependence of the THz emission power 
on the carrier concentration. Furthermore, the detailed 
characteristics of the THz emission confirmed the pres- 
ence of a Dirac-fermion-normalized phonon branch and 
valuable information regarding the fundamental proper- 
ties of Dirac fermions. 
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